Can local bulk effects explain the galactic dark matter? 
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Abstract 



We obtain the virial theorem within the context of a brane-world model without mirror sym- 
metry or any form of junction condition. Taking a constant curvature bulk (neglecting non-local 
00 ■ bulk effects), the local bulk effects generate a geometrical mass, contributing to the gravitational 

energy which may be used to explain the virial mass discrepancy in clusters of galaxies. We fix 
the parameter of this model in agreement with observational data. 
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1 Introduction 

The question of dark matter is one of the most intriguing aspects of the present-day observational 
astrophysics and cosmology. Today, its existence at the galactic and extra-galactic scales is well 
established and is based on two important observations, namely the mass discrepancy in clusters of 
galaxies and the behavior of the galactic rotation curves [1] . Galaxy clusters are the largest virialized 
structures in the universe, and their mass content is supposed to be representative of the universe as 
a whole. The total mass of galaxy clusters ranges from 1O 13 M for groups up to a few 10 15 Mq for 
very rich systems. Cluster masses can be determined in a variety of ways. The application of the 
virial theorem to positions and velocities of cluster member galaxies is the oldest method of cluster 
mass determination [2]. More recent methods are based on the dynamical analysis of hot X-ray 
emitting gas [3] and on the gravitational lensing of background galaxies [4]. The mass determined 
from such dynamical means is always found to be in excess of that which can be attributed to the 
visible matter. This is known as the "missing mass problem." The existence of dark matter was 
not firmly established until the measurement of the rotational velocity of stars and gases orbiting 
at a distance r from the galactic center could be done with reasonable accuracy. Observations show 
that the rotational velocity increases near the center of a galaxy and approaches a nearly constant 
value with increasing distance from the center. The discrepancy between the observed rotational 
velocity curves and the theoretical prediction from Newtonian mechanics is known as the "galactic 
rotation curves problem." This discrepancy is explained by postulating that every galaxy and cluster 
of galaxies is embedded in a halo made up of some dark (invisible) matter. 

The idea that our familiar 4-dimensional (4L>) space-time is a hypersurface (brane) embedded in 
a 5D bulk [5, 6, 7] has been experiencing a phenomenal interest during the last decade. According 
to this brane-world scenario, all matter and gauge interactions reside on the brane, while gravity 
can propagate in the whole 5D space-time. Several brane-world cosmologies have been proposed in 
the context of the Randall-Sundrum (RS) formulations [6], defined in a 5D anti-de Sitter space-time 
(AdSs). The dynamics of these models feature boundary terms in the action and sometimes mirror 
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symmetry, such that bulk gravitational waves interfere with the brane-world motion. This usually 
comes together with junction conditions producing an algebraic relationship between the extrinsic 
curvature and the confined matter [8, 9]. The consequence is that the Friedman equation acquires 
an additional term which is proportional to the square of energy density of the confined matter field 
[10, 11]. This term was initially considered as a possible solution to the accelerated expansion of 
the universe. However, soon it was realized to be incompatible with the big bang nucleosynthesis, 
requiring additional fixes [11]. 

Brane-world scenarios under more general conditions and still compatible with the brane-world 
program have also been rather extensively studied over the past decade where it has been shown that it 
is possible to find a richer set of cosmological solutions in accordance with the current observations [12]. 
Under these conditions, without using Z2 symmetry or without postulating any junction condition, 
Friedman equation is modified by a geometrical term which is defined in terms of the extrinsic 
curvature, leading to a geometrical interpretation for dark energy [13]. There has also been arguments 
concerning the uniqueness of the junction conditions. Indeed, other forms of junction conditions 
exist, so that different conditions may lead to different physical results [9]. Furthermore, these 
conditions cannot be used when more than one non-compact extra dimension is involved. Against 
this background, an interesting higher-dimensional model was introduced in [14] where particles are 
trapped on a 4D hypersurface by the action of a confining potential. The dynamics of test particles 
confined to a brane by the action of such potential at the classical and quantum levels were studied 
in [15]. In [16], the same brane-world model was studied, offering a geometrical explanation for the 
accelerated expansion of the universe. A geometrical explanation for the generalized Chaplygin gas 
was considered in [17] along the same lines. We have also studied exact solutions of the vacuum field 
equations on the brane for two interesting cases. The first solution can be used to explain the galaxy 
rotation curves without assuming the existence of dark matter and without having to resort to the 
Modified Newtonian Dynamics (MOND) [18, 19], and the second solution represents a black hole in 
an asymptotically de Sitter space-time [20]. 

It has been argued that modified theories of gravity based on the RS brane-world scenario can 
explain the observed galactic rotation curves without introducing any additional hypothesis. In 
such models, a spherically symmetric brane-world metric results in a Schwarzschild mass parameter 
receiving a new contribution due to the Weyl tensor which is interpreted as the mass of dark matter 
[21]. Neglecting the cosmological constant, the dark mass increases linearly with the radial distance 
from the galactic center, explaining the flatness of the galactic rotation curves. The bending angle 
of light predicted by brane-world modeles was found to be much larger compared to the predictions 
of dark matter models, the deviation increasing with the degree of compactness [22]. Using the 
smallness of the rotational velocity, a perturbation scheme for reconstructing the metric in a galactic 
halo in brane-world models with induced gravity was developed in [23]. It has also been shown that 
the non-standard model bulk fields can replace dark matter in explaining structure formation, as the 
evolution on the brane becomes similar to that of the cold dark matter [24]. Similar interpretations 
of dark matter as bulk effects have also been considered in [25, 26]. 

As was mentioned above, the galactic dynamics of massive test particles may be understood 
without resorting to dark matter in the framework of f(R) modified theories of gravity [27]. Recently, 
it has been argued that the virial theorem mass discrepancy in clusters of galaxies can be accounted 
for within the context of the RS model and f(R) modified theories of gravity [28, 29]. It would 
therefore be of interest to study the latter in a brane-world model with a constant curvature bulk 
without using the Z2 symmetry or without postulating any form of junction conditions [12, 13]. In so 
doing the gravitational field equations on the brane are modified by a local extra term, Q^ v . Using 
the collisionless Boltzmann equation, we obtain the virial theorem which is modified by an extra term 
that may be used to explain the virial mass discrepancy in clusters of galaxies. 
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2 Field equations for the brane-world 



The embedding of the brane-world in the bulk plays an essential role in the covariant formulation of 
the brane-world gravity, since it tells us how the Einstein-Hilbert dynamics of the bulk is transferred to 
the brane-world. However, there are many different ways to embed a manifold into another, classified 
as local, global, isometric, conformal, rigid, deformable, analytic or differentiable. The choice of one 
or other depends on what the embedded manifold is supposed to do. 

Generally, there are three basic postulates in the geometrical approach considered in brane-world 
scenarios, that is, the confinement of the standard gauge interactions to the brane, the existence 
of quantum gravity in the bulk and finally, the embedding of the brane-world. All other model 
dependent properties such as warped metric, mirror symmetries, radion or extra scalar fields, fine 
tuning parameters like the tension of the brane and the choice of a junction condition are left out as 
much as possible in our calculations [12]. 

In the following we extend this embedding to any constant curvature space, including anti-de 
Sitter AdS$, de Sitter dS$, and the flat M5 cases, with the bulk Riemann tensor written as 

TLabcd = K {QacQbd - QadQbc) , (1) 

where Gab is the bulk metric. In the normal Gaussian frame defined by the embedded space-time, it 
may be decomposed as 

""-("ST L)- "-'-^ (2) 

where k* denotes the bulk constant curvature. In the flat case fc* = and we may write k* = e^g— , 
where is the bulk cosmological constant. We note that for a dS$ bulk e = +1 and for a AdS$ 
bulk e = — 1. With this assumption the Gauss-Codazzi equations [30] 

Ra/3yS = -{K ai K(is - K a sKp 7 ) + TIabcdZ^Z^Z^Z® , (3) 

K a\p-,i\ =T^abcdZ^N b Z'jjZ^, (4) 
which are the integrability conditions for the embedding equations, reduce to 

Ra/3jS = -(Ka~/K(38 ~ Ka5 K fa) + k*(9a-y9/3S ~ 9aS9f3-y), (5) 



Ka\P;i[ = 0, (6) 

where K^ u denotes the extrinsic curvature. From the contractions of equation (5) we can obtain 
the Ricci scalar R of the brane metric, in terms of the Ricci scalar of the bulk 1Z. After adding a 
confined matter Lagrangian and a four dimensional cosmological constant A, the effective Lagrangian 
compatible with the embedding reads, for details see [31] 

C eff = RV=9 - \ (K^K^ + K 2 )^g- + K^g- - - £ ^^9 + - An, (7) 

where K = g^ v K^ u is the mean curvature of the brane-world, y denotes the fifth coordinate and 
1Z = — 20/c* as derived from equation (1). The total derivative term with respect to y can be eliminated 
provided the motion of the brane-world occurs between two fixed minimal boundary hypersurfaces 
where K = [12]. Similar to the RSII model, such boundaries can be moved away so that all 



3 



boundary generated bulk waves are eliminated. Variation of the action with respect to g^ v gives the 
dynamical equation compatible with the embedding and with the confined matter hypotheses 



G flv = %'KGT llv — \g llv + Q llv , (8) 

where 

1 



Quu — 

e 



K^K pu - KKp, - \ (K af3 K^ - K 2 ) 



(9) 



Here, t^ u is the confined matter energy- momentum tensor on the brane and A = — 3fc* + A. As can 
be seen from definition of Q^ u , it is independently a conserved quantity, that is 

= 0. (10) 

so that there is no non-gravitational exchange of energy between this geometrical correction and the 
confined matter. Such an aspect has one important consequence, that is, if Q^ v is to be related 
to dark energy, it does not exchange energy with ordinary matter, much the same as in coupled 
quintessence models [32]. 

In order to solve the Codazzi equation (6), we choose the static spherically symmetric metric on 
the brane in the form 

ds 2 = - e ^ r) dt 2 + e u{r) dr 2 + r 2 (d6 2 + sin 2 6d<p 2 ) . (11) 

Then the York relation 

K - 1 d9lxv C121 

shows that in a diagonal metric, K^ v are diagonal. Here, £ is a small parameter orthogonal to the 
brane that parameterizes the extra noncompact dimension. Now, separating the spatial components, 
the Codazzi equation reduces to 

Kfj,v,a ~ KypT^ = K^ a u — KupT^, (13) 

- (f) *» = - (#) *». <") 

A-J2,i - (i) A' 22 = (re-") A'„. (15) 

The first equation gives ifoo,a = Kn a = K 22>a = #33,(7 = for a = 0,3. Repeating the same 
procedure for a = 2, we obtain Kqq^ = K\\^ = K 22>a = 0. This shows that K\\ depends only on 
the variable r. Assuming K\\ = ae^ r ' and using equations (14) and (15), one finds 

K m {r) = -ae^ + ce^/ 2 , (16) 

K 22 {r) = ar 2 + (3r. (17) 
Taking (i, v = 3 in the first equation we obtain 

#33,1 " (J) #33 = (e~ u r sin 2 o) K u = ar sin 2 9, (18) 
#33,2 - (cot 9)K 33 = (sin 9 cos 6) K 22 . (19) 
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Using equations (16), (17) and (18) then results in 

^33 ( r i 6) = ar ~ 2 sm2 @ + r/3sin 2 8 + rc\ sin 8. 
Now, use of equation (9) leads to the components of Q^ v 

3a 2 r 2 + 4a(3r + (3 2 + -^(2ar + /?) 



(20) 



Qoo = 


9oo 






£r 2 


Qu = 


9u 




er 2 


Q22 = 


922 






er 


Q33 = 


933 - 


er - 



sin 8 

3a 2 r 2 + 4a/3r + (5 2 + f 2ar + /? - ere"' 1 / 2 ) - 2ce~^ 2 (ar 2 + (3r) 
sm.8 \ / V / 



sin# 



(21) 



Since G\ = G| and thus Q\ = Q% one obtains C\ = 0. Now, using these relations and equation (8), 
the gravitational field equations become 



^ "^2+7 +^ = 8 ^// + A + -(3a 2 + 



(22) 



e"" f 4 + ±\ ~ ~ 2 = **G£) f - A + - (-3a 2 + 2ace-/ 2 - ^ - ^ + , (23) 

«- - £ - *f + / + f ) = ™>ti) - » + I (^ + ^»-" /2 - ^ + T 2 '-'" 2 ) .('24) 

where a prime represents differentiation with respect to r and e is the signatures of the bulk space, 
e = +1 for the (IS5 bulk and e = — 1 for the AdS§ bulk. In the above equations we have assumed that 
the matter on the brane consists of an anisotropic fluid, characterized by an effective energy density 
p e ff / 0, a radial pressure p^jj and a tangential pressure Pg//- 

In the next section, we will investigate the influence of the extra term, Q^u, on the dynamics of 
the galaxies by adopting the method introduced in [28] . 



3 Dark matter and the role of extrinsic curvature 

In order to obtain the virial theorem for galaxy clusters in our model, we have to write down the gen- 
eral relativistic Boltzmann equation governing the evolution of the distribution function /. Galaxies, 
which are treated as identical and collisionless point particles, are described by this distribution func- 
tion. For the static spherically symmetric metric given by equation (11) we introduce the following 
frame of orthonormal vectors [33, 34] 

ef = e^5% e W = e "% ef = r6 2 , ef = rsm8S% (25) 

where g^efflefi = r/ a )( b ). The four velocity of a typical galaxy, satisfying the condition v^v^ = 
— 1, in tetrad components is given by 

=v ^e^\ (26) 
The relativistic Boltzmann equation in tetrad components is given by 

M e P ^L + J a ) Jb)Jc) d f _ Q (27) 
( a W + 7 ( 6 )( c ) dvW ~ ' { ' 
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— e^le^,,e?„s are the Ricci 



where the distribution function / = f(x^,v^), a = 0,1,2,3, and 7( 6 )( c ) — c p;o- c (f,) c ( c ) 
rotation coefficients. Thus, equation (27) becomes 



Or 



' v\_dp_ 

r 86 ' rsin# dp \ 2 dr 



e v l 2 v 6 df e u / 2 v v df 



df_ 

dv r 



Vr df Of 

— I v e^r- + v f-- 



W2» 



df 



df 



dv e 



dv. 



cot 6 v$ — Vp— — 



dv, 



dv e 



0. 



(28) 



Assuming that the distribution function is only a function of the radial coordinate r, equation (28) 
reduces to 



df ( v\ dn j+jT 



dr 



2 dr 



df_ _Vr_ ( df_ +v df_ 

dv r r \ dve v dv v 



W2, 



df 



— cot ( Vg V 



dvu 



df 

^ dv e 



0. (29) 



Now, taking into account the spherically symmetric nature of the problem at hand and integrating 
over the velocity space and assuming that the distribution function vanishes rapidly as the velocities 
tend to ±oo, one obtains 



d_ 

dr 



P(v 2 r ) 



1 

+ 2 P 



(v 2 ) + (v 2 ) 



dfi 
or 



(vl) + (v 2 ) - 2(v 2 ) 



0. 



Multiplying equation (30) by 4irr 2 and integrating over the cluster leads to 
rR 



1 f R r 
r 2 dr + - / 4nr 3 p (v 2 ) + (v 2 ) 

2 Jo L 



-jT 47TP [(V 2 ) + (V 2 ) + (Vl) 

Now, using the total kinetic energy of the galaxies 

k = j\i:p[(v 2 .) + (v 2 e) + (vl) 

equation (31) is reduced to 

2k- l -j\^p[(v 2 ) + (v 2 r ) 



—dr = 0. 

Or 



r 2 dr, 



%r = 0. 

Or 



(30) 



(31) 



(32) 



(33) 



Let us write the energy-momentum tensor of the confined matter on the brane in terms of the 
distribution function 



1 jlV 



= J fmv^v u dv, 
which gives 

Peff = P(Vt), Peff = P( V rX P e )) = p( v l) = pltfp) ■ 

Now, using these relations and summing the gravitational field equations (22)- (24) we find 



-v I P_ _ , P_ , M 
3 r 4 + 2 + 4 



'pj = 4nGp(v 2 ) -X + - £ (-3a 2 - ^ + Sace^ 2 + ^e"^ 2 ) 



(34) 



(35) 



(36) 



where (v 2 ) = (v 2 ) + (v 2 ) + (i>|) + (v 2 ). In order to obtain the generalized virial theorem in our model 
we have to use some approximations and assumptions. Since the dispersion of the velocity of galaxies 
in the clusters is of the order 600 - lOOOkm/s, that is (f ) 2 «4x 10~ 6 - 1.11 x 10~ 5 <C 1, therefore 
we can neglect the relativistic effects in the Boltzmann equation, and use the small velocity limit 
approximation, so that (v 2 ) ~ (v 2 ) ~ (v 2 ) <C (v 2 ) ~ 1. The intensity of the gravitational effects can 
be estimated from the ratio of the Schwarzschild radius to the radius of cluster, GM/R, which for 
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typical clusters is of the order of 10" 6 <C 1. Therefore the gravitational field inside galactic clusters 
is weak and we can use the weak gravitational field approximation, so that the quadratic terms can 
be neglected [28]. Thus, equation (36) is given by 

1 d (_^d\t 

where 



4nG Pextr = - (-3a 2 - ^ + Sace^' 2 + ^e^ 2 ) , (38) 
e \ r r J 

and the suffix 'extr' stands for 'extrinsic'. Multiplying equation (37) by r 2 and integrating from to 
r we obtain 



1 2^ 

2 r dr 



,-_ _ GM(r) + l\r 3 - 4vrG f p extr {r')r' 2 dr' = 0, (39) 
3 Jo 



where M(r) is the mass out to radius r, so that dM(r) = \-npr 2 dr. The total baryonic mass of the 
system is given by M = JJ* Airpr 2 dr'. We also define the geometrical mass as 

M extr {r) = 4tt / p extr {r')r' 2 dr'. (40) 
J o 

Multiplying equation (39) by dM(r) and integrating from to R and introducing the moment of 
inertia of the system as I = J R r 2 dM(r), we obtain 

2k + W + ^XI - W extr = 0. (41) 

In the above equation we have used equation (32) and the following definitions 

r R GM(r 



W -- 
and 

W ex tr 



-f^MH, (42) 
Jo r 

Jo r 



where W is the gravitational potential energy of the system. For the case a = (3 = c = 0, we obtain 

W ex tr = and equation (41) reduces to the virial theorem in the presence of a cosmological constant 
[33]. Now, let us introduce the radii R v , Rj and R ex tr as [28] 

M 2 

Rv = r R M(r)dM(r) i ( 44 ) 
JO r 

9 \*r 2 dM(r) , , 
M 2 xtr 

Rextr = f-R M extr {r)dM(r) ( 46 ) 
JO r 

where R v is the virial radius and R ex tr is defined as the geometrical radius of the clusters of galaxies. 
Defining the virial mass as [33] 

a-^, (47) 
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and using the following relations 



W = W extr = ^^, I = MRl (48) 



the generalized virial theorem (41) is simplified as 



1/2 



M 



3GM ^ M2R extr ) ' 



As can be seen, we have three type of mass in equation (49), namely, the total baryonic mass of the 
system represented by M (including the baryonic mass of the intra-cluster gas and of the stars, other 
particles like massive neutrinos), the virial mass represented by M v and finally, the geometrical mass 
represented by M extr . 

On large distance scales associated with galaxies, we can ignore the contribution of the cosmo- 
logical constant relative to the mass energy of the galaxy Since, as it turns out, M v is considerably 
greater than M for most of the clusters, we can neglect the unitary term in equation (49). Therefore, 
the virial mass is given by 



M v ~ M ex tr 



R X 1 / 2 

-) • (50) 



Rextr J 



This equation shows that the virial mass is proportional to the geometrical mass due to the local extra 
term, Qui/, thus the virial mass discrepancy in clusters of galaxies can be explained by the generalized 
virial theorem. Since galactic clusters are dark matter dominated objects, the main contribution to 
their mass comes from the geometrical mass M extr so that with a very good approximation M extr ~ 

Mtotal ■ 

All the way, among various possible solutions that depend on different choices of arbitrary con- 
stants a, (3, c, we consider a = and (3, c ^ 0. Since the bulk of matter in clusters of galaxies is in the 
form of dark matter, we can neglect the effect of the ordinary matter in equations (22)- (24). Thus, 
the solutions are given by 



e- = 1 - 0- , (51) 




e -,/2 = Oll« I. (52) 

pc r 



Now, using equations (38) and (40) the geometrical mass is given by 

r 2 dr . (53) 



GM extr {r) = - f 
£ Jo 

Substituting e~^/ 2 from equation (52), we obtain 



3a 2 -^ + 3ac e -^')/ 2 + ^e-^')/ 2 



GM extr (r) = 2 ^1 - r. (54) 

Since according to our physical interpretation M extr is an effective geometrical mass it must satisfy 
the condition M extr > 0, therefore we take e = — 1 (AdSs bulk). This geometrical mass is linearly 
increasing with r, thus having a similar behavior as that of dark matter in clusters of galaxies. 
An important observational quantity is the radial velocity dispersion a 2 , which is related to the 
total mass in a cluster by the relation GM V (R V ) = a^R v . For r = R v equation (54) becomes 
GM extr (R v ) ~ (1 + (3 2 )Rv and assuming that M extr ~ M v we obtain 

(l + /3 2 )«a r 2 . (55) 
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Therefore, the present model gives a geometrical interpretation for the mass discrepancy in clusters 
of galaxies, supporting our previous results that the galaxy rotation curves can be explained without 
assuming the existence of dark matter and new modified theories (Modified Newtonian Dynamics) 
[20]. 

Finally, we should mention here that there is a difference between our model and the work pre- 
sented in [28]. In this paper we have shown that the virial mass of galaxy clusters is mainly determined 
by the geometrical mass associated with the extrinsic curvature terms (local bulk effects) whereas in 
[28], the virial mass is specified by the dark mass associated with the transmitted projection of the 
bulk Weyl tensor (non-local bulk effects). 

4 Conclusions 

In the RS brane-world models, the 4-dimensional effective Einstein equation has some extra terms 
known as dark radiation and dark pressure, which describe the non-local bulk effects due to the 
gravitational field of the bulk. It has been shown that the dark radiation term can be used to explain 
the virial theorem mass discrepancy in clusters of galaxies [28] . 

In the present paper we have generalized the virial theorem within the context of brane-world 
models without using the Z2 symmetry or without postulating any junction condition. To obtain the 
virial theorem we have used a method based on the collisionless Boltzmann equation. Assuming a 
constant curvature bulk, the virial theorem was modified by a local extra term, equation (41), which 
was then used to explain the virial mass discrepancy in clusters of galaxies. Thus, the geometrical 
mass can play the role of dark matter in clusters of galaxies, supporting our previous results that 
the rotational galactic curves can be explained in our model without introducing any additional 
hypothesis [20]. Finally, after fixing the parameter of this model, we found them to be in good 
agreement with observations. 
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